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Abstract. We discuss factorization in heavy quarkonium production in high energy collisions using 
NRQCD. Infrared divergences at NNLO are not matched by conventional NRQCD matrix elements. 
However, we show that gauge invariance and factorization require that conventional NRQCD 
production matrix elements be modified to include Wilson lines or non-abelian gauge links. With 
this modification NRQCD factorization for heavy quarkonium production is restored at NNLO. 


INTRODUCTION 


Quarkonium production and decay have been the subject of a vast theoretical literature 
and of intensive e^merimental study (ll], in whieh the effeetive field theory nonrelativistie 
QCD (NRQCD) [Q] has played a guiding role. NRQCD offers a systematie formalism 
to separate dynamies at the perturbative mass seale of the heavy quarks from nonpertur- 
bative dynamics, through an expansion in relative velocity within the pair forming the 
bound state. An early sueeess of NRQCD was to provide a framework for the Tevatron 
Run I data on high-pr heavy quarkonium produetion fsl], and it has been extensively 
applied to heavy quarkonia in both eollider and fixed target experiments HI] . 

In eontrast to quarkonium decay, fully convincing arguments have not yet been given 
for NRQCD faetorization as it is applied to high-py produetion processes. Here we 
review new tests of NRQCD faetorization for high-pr quarkonia produetion at next-to- 
next-to-leading order (NNLO) 10] • At NNLO we find infrared divergences that do not 
fall precisely into the pattern suggested in Ref. 10] • These divergences may, however, be 
ineorporated into eolor oetet matrix elements by a teehnieal redefinition in whieh Wilson 
lines whieh makes them gauge invariant, restore faetorization at NNLO. 

If we assume NRQCD faetorization, we have 
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where the are NRQCD operators for state H, whieh were introdueed in Ref. |0] in 
the form 
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FIGURE 1. Representative NNLO contributions to g ^ cc fragmentation in eikonal approximation. 


where ah is the ereation operator for state H,xiW) two eomponent Dirae spinors and 
where jtn and involve produets of eolor and spin matriees, and at higher dimensions 
of eovariant derivatives. 

Our gauge invariant redefinition of produetion operators in oetet representation is the 
replaeement (whieh we refer to as a gauge eompletion), 

(0) ^ Z^X,cV^(0)^?[0,A]ez, [ahan) , (3) 

where we have exhibited the eolor indiees of the oetets. Here 

noo 

4>/[0,A] = ^exp[—/g / dXl-A{Xl)] (4) 

Jo 

is a Wilson line or gauge link, with 1^=0. The field A=A^T^ with being the generator 
in the adjoint representation. 

Representative diagrams for the fragmentation funetion at NNLO are shown in Fig. 1 
in whieh the off shell lines are indieated by the heavy dots. The full infrared fragmen¬ 
tation funetion is generated by taking all allowed euts of the remaining lines of eaeh 
sueh diagram. We are eoneerned only with diagrams that eonneet oetet to singlet quark 
states, and whieh are not topologieally faetorized, sinee these are the potential sourees of 
nonfaetoring behavior in both the fragmentation funetion and related eross seetions. The 
original argument for NRQCD faetorization was based on the eonjeeture that all infrared 
regions in these diagrams eaneel after this limited sum over euts ||2]- In faet, this is the 
ease at NNLO only if we employ the gauge-eompleted definitions for NRQCD matrix 
elements, as in Eq. © above. 

The individual elasses of diagrams in Fig. la and lb, for whieh two gluons are 
exehanged between the quarks and the Wilson line, satisfy the infrared eaneellation 
eonjeeture of Ref. 10], by summing over the possible euts and eonneetions to quark 
and antiquark lines, as do diagrams that have three gluon-eikonal vertiees on the quark 
pair and one on the Wilson line. For the elass of diagrams related to Fig. le, however, 
with a three-gluon interaetion, this eaneellation fails. Expanding to seeond order in the 
relative momentum q, the full eontribution from Fig. le, found by eutting the gluon line 
k\ and the Wilson line, ean be written as 
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where Vy^x [ki, ^ 2 ] represents the momentum part of the three-gluon vertex. Summing 
over all eontributions, however, we find a noneaneelling real infrared pole in the frag¬ 
mentation funetion (in the quarkonium rest frame P =0) 0], 
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where v is the relative veloeity of the heavy quark pair. 

Eq. ® shows explieitly the breakdown of topologieal faetorization of infrared de- 
pendenee at NNLO. Its presenee implies that infrared poles would appear in eoeffieient 
funetions at NNLO and beyond when the faetorization is earned out with oetet NRQCD 
matrix elements defined in the eonventional manner, Eq. Q. On the other hand, when 
defined aeeording to its gauge-eompleted form ©, eaeh oetet NRQCD matrix element 
itself generates preeisely the same pole terms given in ® above. Thus, at least at NNLO 
and to order v^, when the non-perturbative NRQCD matrix are defined aeeording to Eq. 
©, the faetorization is restored. We are investigating the extensions of this result. 
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